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I-  HISTOfil  AND  INTRODOCTION 


This  thesis  is  the  third  in  a  series  of  investigations 
into  a  proposal  put  forth  by  E.  Lauer  of  Naval  Underwater 
Systems  Center,  New  London,  CT,  in  1979.  In  his  memorandum 
£Bef.  1],  Lauer  describes  a  "new”  way  to  analyze  sound 
propagation  in  the  ocean  and  two  applications  of  the  method, 
which  he  named  the  "Navenuaber  Technique,"  or  F. T-  In  his 
proposal,  Lauer  stated  that  a  source  of  sound  in  the  ocean 
could  be  pinpointed  in  both  depth  and  range  by  a  single 
omnidirectional  hydrophone  provided  the  source  is  generating 
a  continuous  wave  tone. 

The  N-T.  was  initially  described  by  F.  DiNapoli  [Bef.  2] 
as  an  intermediary  step  in  his  development  of  a  speedier 
computer  algorithm  used  in  analyzing  sound  propagation  as  a 
function  of  range.  CiNapoli  and  Lauer  proposed  converting 
the  pressure  as  a  function  of  range  ^(r)  to  the  pressure 
spectrum  as  a  function  of  the  wavenumber  P(k).  The  conver¬ 
sion  is  accomplished  by  taking  a  weighted  Fourier  transform 
of  the  pressure  field  ^(r).  Once  the  pressure  spectrum  is 
obtained,  analysis  of  the  pressure  field  can  be  accomplished 
in  a  manner  analogous  to  that  presently  used  in  signal 
processing. 

Lauer  proposed  two  uses  for  his  "Savenumber  Technique:" 

1.  determination  of  the  depth  of  the  acoustic  source,  an 
ability  which  has  obvious  tactical  applications; 

2.  use  of  the  K.  I.  to  evaluate  the  accuracy  of  existing, 

and  future,  acoustic  models  such  as  P.E.  and 

F.A.C.T.,  by  breaking  the  generated  pressure  field 
down  into  ray  path  "families,"  such  as  bottom-bounce, 
ref racted-surface-ref lected,  and  surf ace- ducted  fami¬ 
lies,  thus  giving  a  quantitative  read-out  of  what 


proportion  of  the  total  acoustic  energy  is  being 
channelled  through  the  various  ray  paths. 

For  both  uses,  a  knowledge  of  the  acoustic  environment  is 
necessary. 

It  is  the  intended  purpose  of  this  thesis  to  investigate 
the  validity  of  the  first  application  of  Lauer's  wavenumber 
technigue  by  using  it  in  conjunction  with  an  acoustic 
phenomenon  that  is  well-known  and  for  which  accurate  results 
can  be  calculated  with  precision;  this  phenomenon  is  the 
LLoyds  Hirror  for  an  Acoustic  Doublet.  Osing  several 
source/receiver  combinations,  this  investigator  intends  to 
compare  actual  source  depth  with  that  predicted  by  the  a. I. 
If  Lauer’s  technigue  is  cogent,  the  source  depths  calculated 
via  the  two  methods  should  be  equal,  or  very  close  within  a 
statistically  acceptable  degree. 

As  mentioned  in  the  beginning  of  this  thesis,  two  prior 
investigators,  B.  Stamey  and  J.  Blanchard,  looked  into 
applying  the  W.  T.  to  determining  source  depth.  Staley's 
investigation  [Ref.  3]  utilized  a  parabolic  equation 
computer  model  developed  by  H. K.  Brock,  the  Split  Step  Fast 
Fourier  Transform,  or  SSFFT,  to  generate  the  pressure 
fields.  He  qualitatively  compared  the  model  to  a  Normal 
Mode  model  and  a  P.  S.  Finite  Difference  model,  both  range 
dependent,  and  to  DiNapoli's  Fast  Field  Program  model,  a 
range  independent  model.  Stamey  concluded  that  the  W.T. 
showed  promise  as  an  analysis  tool  and  that  further  investi¬ 
gations  were  warranted. 

J.  Blanchard  carried  Stamey's  researches  one  step 
further  [Ref.  4].  Using  two  Parabolic  Equation  computer 
algorithms.  Brock’s  Split  Step  Fast  Fourier  Transform  and 
Jaeger’s  Implicit  Finite  Difference,  as  his  pressure  field 
generators,  he  examined  the  use  of  the  spacing  between  nulls 
of  the  pressure  spectrum  to  determine  source  depth.  His 
results  are  interesting  and  support  the  need  for  further 
invest  igation. 


The  author  was  not  satisfied  with  Stanley's  and 
Blanchard's  findings.  Both  individuals  acknowledged  short¬ 
comings  in  their  respective  studies,  especially  with  regard 
to  the  mutual  presence  of  the  "0-shaped  phenomenon"  encoun¬ 
tered  in  the  pressure  spectrum.  One  should  remember, 
however,  that  their  theses  were  preliminary  studies  only  and 
were  produced  within  a  highly  restricted  tine  frame.  Also, 
they  used  existing  computer  algorithms  specifically  designed 
to  approximate  acoustic  propagation  in  a  velocity-variant 
medium.  However,  the  environmental  settings  used  by  both 
gentlemen  were  that  of  the  Lloyds  Mirror  for  an  Acoustic 
Doublet  which  requires  a  constant  sound  speed  throughout  the 
water  column;  thus,  the  approximations,  assumptions  aid 
"fudge  factors"  used  by  these  models  make  their  results  for 
an  isospeed  medium  highly  artificial  and  somewhat  suspect. 
Fortunately,  however,  there  exists  for  the  Lloyds  Mirror  a 
simple,  geometric  solution  specifically  designed  for  an 
isospeed  environment  [Ref.  5].  Since  the  pressure  fields 
for  varying  source  and  receiver  depths  and,  therefore,  the 
corresponding  pressure  spectra,  can  be  precisely  calculated, 
it  was  thought  this  model  would  be  a  good  check  on  the  oper¬ 
ational  applicability  of  Lauer's  proposed  technique. 

Secondary  considerations  of  this  study  were  to  examine 
what  effect,  if  any,  the  introductions  of,  firstly,  surface 
waves  and,  secondly,  noise  would  have  on  the  pressure  spec¬ 
trum.  The  first  objective  »'as  simulated  by  use  of 
mathematical  formulae  given  in  reference  5;  for  the  second 
objective,  it  was  thought  that  complex  noise,  random  in  both 
amplitude  and  phase  and  similar  to  a  Rician  distribution 
[Ref.  6;  p.  189],  would  provide  a  simple  but  reasonable 
approximation  of  an  ocean  noise  field.  This  investigator 
intends  to  gradually  intensify  the  surface  waves  and  the 
noise  field,  separately,  until  the  original  pressure  spec¬ 
trum  is  no  longer  recognizable.  In  this  manner,  it  should 


be  possible  to  make  a  qualitative  assessment  of  their 
respective  debilitating  effects  under  the  carefully' 
controlled  conditions  found  in  the  Lloyds  Mirror  phenomenon. 

Section  II  sets  forth  in  more  detail  the  theoretical 
development  of  each  point  described  in  this  introduction. 
Section  III  sunimari2es  the  investigator's  results  and  the 
conclusions  drawn  frcm  those  results.  Section  IV  contains  a 
listing  of  the  computer  algorithms  utilized. 
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II.  THEORI 

This  chapter  presents  the  mathematical  basis  of  the 
wavenumher  technique. 


A.  TBE  LLOIDS  HIEBOB  PHEHOHENON 

A  diverging  mcnoxrequency  spherical  pressure  wave 
[Ref-  5:  p.112]  can  be  written  in  complex  form  as 


-//i#  A  ^i(kR-(x}f) 
£(R,f)-£<R>e-"^'  =  -pS 


(egn  2.1) 


where  ^(E)  is  the  spatial  factor.  Referring  to  Figure  2.1, 
Hj  is  the  range  from  the  source  to  the  receiver  and  32,  is 
the  range  from  the  image  to  the  receiver.  For  convenience 
in  all  tha^  follows,  we  shall  set  A  to  unit  magnitude.  The 
following  equations  apply  to  the  acoustic  waves  propagating 
directly  out  from  the  source  JPj  and  from  the  apparent  image 
P-  (which  is  actually  the  wave  reflected  from  the  surface) : 


=  4(R)e 


(egn  2.2) 


where  P  #  P,  ,  and  P?  are  complex  functions  of  horizontal 
range  r, vertical  depth  Z  and  time  t;  the  minus  sign  in  the 
equation  for  P;  is  derived  from  the  surface  reflection  coef- 
ficient,  =  -1,  for  a  smooth  surface.  Inspection  of 
Figure  2.1  reveals  that 


18 


Rt.  PRODUCED  AT  GO  V  r  r?  rVf  t 


id  hr  =  R-R, 


And  so  the  total  field  can  be  written  as 


-P(R)c‘^'-  ^-5 - 

'v 


i{kRr  (jJfJ 


(eqn  2.3) 


Since  P.  and  P;  both  have  the  same  time  factor, 
exp  (— (.wt) ,  we  can  retain  just  the  spatial  factors  and  equa¬ 


tions  2.3  reduce  to 


PfR)  = 


r  _ 

L  p, 


e‘^^2 

Ps. 


(eqn  2.4) 


Zqudtion  2.4  is  a  form  of  the  complex  pressure  as  a  function 
of  range  used  in  the  computer  algorithm  shown  in  Appendix  I. 

Inspection  of  Figure  2.  1  reveals  that,  for  R  >>  2^  and 
Q  very  small, 

r^R 


sin  &  ^ 


(eqn  2.5) 


aiid  the  pressure  amplitude  can  be  approximated  by 


2  fkZ.Z^ 


(eqn  2.6) 


Looking  at  just  the  formula  for  the  pressure  amplitude, 
we  can  see  that  as 


REPRODUCED  AT  GOVrT’r'if’rtJf 


nrc  i  n-=.0jlj2y3j.. . 


(eqn  2.7) 


the  pressure  amplitude  goes  to  zero,  producing  the  classic 
|P(E)  1  VS  E  curve  shown  in  Figure  2.2 


B.  THE  EELATIOHSHIP  EETSEEN  K,  y  ,  AMD  Q 

In  his  description  of  the  Ravenumber  Technigue  [Eef.  1: 
p.  5-6],  Lauer  utilizes  the  horizontal  and  vertical  compo¬ 
nents  of  the  wavenumber,  k.  The  general  relationship  amoung 
these  three  terms  is  illustrated  in  Figure  2.3  and  can  be 
written  mathematically  as 

k  =  the  wavenumber  *  2vf/c 

'y  =  the  horizontal  component  of  k  *  k  cos(p 
P  =  the  vertical  component  of  k  *  k  sin  Cp 
or 


.  k  -\/r^ 


2 


(egn  2,3) 


C.  THE  RAVERaaBEE  TECHHIQOE  AND  THE  LLOYDS  HIESOE 

I  Given  a  point  source  in  free  space,  the  monof reg uency 

pressure  field  at  the  receiver,  ^(H)  ,  can  be  expressed  as  a 
spherical  wave  (time  factor  suppressed). 


I 
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£(R}  =  i 


ikR 


(eqn  2.9) 


where,  from  Figure  2.1, 


,  for  >  Z^ 


or 


=  r^  +  Z  ,  for  Z  <  Z 
s  r  n 


(eqn  2.  10) 


In  integral  form,  ^(E)  can  be  written  as  [fief.  7:  p.  127] 

“  J>(7r) 


P(R)  =  1 


ikR 


r< 

Jr\ 


iP 


7  dy 


(egn  2.11) 


which  is  taken  from  the  Fourier-Bessel  Transform  Fair 
[Ref.  7:  p.  126], 

l(r)  =  j'^g(7)Jj7r)7d7 


^(y)  =  f{r)Jj7r)rdr  2.  12) 

where  represents  the  acoustic  pressure  function  in  the 
range  domain,  and  ^('y)  is  the  acoustic  pressure  spectrum  in 
the  wavenumber  domain.  The  sign  of  the  exponential  function 
in  the  integral  is  based  on  which  of  Zj  and  is  greater. 
For  the  case  of  receiver  depth  being  greater  than  source 
depth  so  that  waves  from  both  image  and  source  are  travel¬ 
ling  downward  at  the  receiver  depth,  the  total  pressure  at 
the  receiver  can  be  expressed  as 


-  J  iQ  d,(7r)7dy  2.13) 


X  "t  fjsr 


Use  of  Euler’s  Identity  reduces  Equation  2.13  to 
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£(R)  =  -2  r 


^  L^W^P^s.)..  ( yr)  y  ay 


(egn  2.  14) 


Using  the  relationship  between  the  HanJcel  functions  [Bef.  5; 
p.  449], 

H['^(yr)  =  JJyr)  +  iYJyr) 


=  JJyr)  -  YJyr) 


(egn  2.  15) 


one  can  re-write  the  Eessel  function  in  eguation  2.14  as 


JJyr)  =  i[H^J\yr)  ^  H!,%r)] 


(egn  2 . 16) 


so  that 


Zg<y)  “  l’°iMH'."(yr)rdr  *  fJi(rJH!%r)rdr 


Letting  r’  =  -r  and  looking  at  the  second  term  on  the  right 
hand  side  of  eguation  2.  17, 

Jr-OO  fo]  r*  ft) 

M-r')  H,  (-yr')r'clr'  =  j  XMH,  (yr)rdr  (egn  2.13) 


Now,  assuming  that  f  (r)  =  f  ("C)  ,  then 


/CO  //I 

Z(r)  H,  (yr)rdr 

oo 


Therefore, 


(egn  2.  19) 
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-  rrnVj. 

-J-00  R 


(egn  2.  20 ) 


For  yr  >  2  7r  ,  the  asymptotic  approximation  of  th< 
Rankel  function  for  large  argument  can  be  used: 


H!%r)  =» 


yf'  (egn  2.  21) 

assuming  |r^(r)|  goes  to  zero  faster  than  (In  r). 


2 

TT  f 

'/T" 

^  Xoo 

.  P 

terms 

inside 

the 

sinfpzj 


e''^'’dy(eqn  2.22) 


defined  as  ^  then 

£(R}  =  -  7^  e-'^[^0yje‘'^'dr 

and 


_f”grTJe''>"'dT  =  E(R)  2. 


23) 


or,  using  eguation  2.12, 


P<R)e-i'f'dr 


(egn  2.24) 


Note  that  equation  2.24  has  the  form  of  the  Fourier  trans¬ 
form  cf  the  pressure  function  if  we  define 


c  f”"  Ncr. 


fjr)  H  /T~ 


(egn  2.25) 


a  cr  rtvjLJ  Ul.  t  U  U  C’ V  L  f  >  rj  f  '  r  |  j  7  ryprrjc 


SO  thatyr  ?{E)  and  gi'y)  are  Fourier  Transform  pairs. 
This  equation  can  be  easily  evaluated  [fief.  8], 


^(y)  =  - 


r  tti 

p 

/ 

/T* 

e 

V2  7r  1 

sin(Pz^) 


(eqn  2.  26) 


for  ^  magnitude  of  ^('y)  is 


I  <T)| 


=  sin(pZs) 


(eqn  2.27) 


For  the  case  where  >  Zj.  >  it  can  be  shown  taat 

TT 


g{y)=  - 


/  4 


and  the  magnitude  of  ^(y)  is 


iPzs.  .  rn  \ 
P  e  sin{JJz;^ 


(eqn  2.23) 


9(r} 


=  ^^sin(0z,l 


( eqn  2. 29) 


A  close  look  at  the  magnitude  cf  the  pressure  spectrum  func¬ 
tion  in  equation  2.27  reveals  that  nodes  occur  for  values  of 


6Z  =  niT  .  n  =  0.1,2,3.  .  .  . 


(eqn  2.30) 
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If  the  magnitude  of  the  pressure  spectrum  is  plotted  as  a 
function  of  the  vertical  wavenumber  then  the  spacing  between 
nodes  A/3  is  uniform  and  the  depth  of  the  source  can  be 
derived  from  the  relationship  given  by 


(egn  2.  31) 

In  the  event  the  source  is  deeper  than  the  receiver,  looking 
at  equation  2.29,  we  can  see  that  nodes  now  occur  for  values 
of 


SZ  -  riTT 
r 


n  =  0,1, 2. 3,... 


( egn  2.  32) 


Plotting 

umber, 

receiver 


equation  2.29  as  a  function 
the  spacing  between  nodes 
depth  based  on 


of  the  vertical  waven- 
will  now  reveal  the 


Z 


r 


TT 

'Kb 


(egn  2.  33) 


Notice  in  equations  2.31  and  2.33  the  exchange  of  source  and 
receiver  depths.  By  placing  the  receiver  at  a  shallower 
depth  than  the  source,  no  new  information  (namely  the  source 
depth)  is  to  be  found. 

In  summary,  there  is  a  way  to  calculate  the  complex 
pressure  amplitude  as  a  function  of  range  using  the  complex 
pressure  spectrum.  Conversely,  if  the  complex  pressure  ^(r) 
has  been  measured  at  the  receiver,  then  the  complex  pressure 
spectrum  ,P  (k)  can  be  derived  by  making  use  of  the  relation¬ 
ship  between  the  Fourier-Bessel  Transform  Pairs  (see  equa¬ 
tion  2.12),  and  the  depth  of  the  source  can  be  found  from 
the  magnitude  of  that  spectrum  (see  equation  2.3  1). 
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The  foregoing  is  the  aathematics  involved  xn  deriving 
the  theoretical  pressure  spectrum.  One  purpose  of  this 
thesis  was  to  compare  the  theoretical  value  of  the  pressure 
spectrum  as  derived  with  the  Fourier-Bessel  Transform  with 
the  pressure  spectrum  obtained  by  use  of  the  computerized 
Discrete  Fourier  Transform  (DFT)  ,  otherwise  known  as  the 
Fast  Fourier  Transform,  or  FFT.  At  this  point,  it  will  be 
helpful  to  review  what  happens  in  the  FFT. 

Starting  with  the  Fourier-Bessel  transform  pair  shown  in 
euuation  2.12,  f  (r)  and  re-written  in  terms 
of  a  discrete  sum  as  is  done  in  the  DFT: 


2'Krr)n 


I  ~  •  g(y)  =  ^  g_(nl\y)  e  =  ^(m/\r) 

niQ 


(egn  2.  34 ) 


where  N  -  1,2,3,-..,  up  to  some  large  positive  integer,  and 
represents  the  number  of  transform  points,  and 


yr  =  (n^y)  (m  Lr)  = 


(egn  2.  35) 


based  on  the  relationship  between  Ar  and  A*y  as  described 
in  sampling  theory  [  Hef .  2:  p.  2], 

ArAr 


A/'  = 


(egn  2 .  36) 


~  NLr 

By  convention,  the  sign  of  the  exponential  function  in  equa¬ 
tion  2.34  is  taken  as  negative  when  performing  the  "forward" 
transform , 


FFT  l(m^r)  =  g(n^y) 


(egn  2.37) 
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and  is  taken  as  positive  when  performing  the  "inverse" 
transform. 


FFT'' 


£(nl\y) 


(egn  2.38) 


Then,  r<»  r  Hd 

/  g(y)  dy  ^ 

CD 


Y^g(n!ly)e‘  ^ 


L  «=/ 


Ay 


where 


gJnLrf)  =  sin(pzs) 


Therefore 


and 


£(mLr)  =  -  ^ 


hy  FFT 


-/ 


g(n  Ay) 


g(n  Ay) 


yfT 

=  — ^  sin(pZs) 


(egn  2.39) 


THE  EFFECTS  OF  SUBFiCE  BOUGHHESS 


Surface  roughness  [Eef.  5:  p.409]  reduces  the  Lloyds 

Mirror  effect.  The  rougher  the  surface,  the  greater  the 
effect  for  angles  increasingly  closer  to  grazing  incidence. 
If  M  represents  the  surface  roughness  factor,  then  equation 

2.3  can  be  re-written  as 

ikR 


P(R)  = 


m2  ikR 
-M  e 


R 


where 


M 


s  sintji 

A 


(egn  2.  40) 
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is  the  grazing  angle  of  incidence,  H  is  the  average 
height  of  the  surface  gravity  waves,  and  \  is  the 
wavelength  of  the  narrowband  continuous  wave  acoustic 
signal.  When  M  is  less  than  1,  the  surface  is  said  to  be 
smooth;  when  M  is  greater  than  1,  the  surface  is  said  to  be 
rough. 

Looking  at  Figure  2.1,  we  can  see  that  the  overall 
effect  of  increasing  surface  roughness  is  to  reduce  the 
contribution  of  the  surface  reflection  or  "image"  to  the 
interference  pattern  at  the  receiver.  In  other  words,  as  H 
increases,  the  energy  detected  by  the  receiver  is  increas¬ 
ingly  only  the  energy  coming  directly  from  the  source.  When 
the  effect  of  surface  roughness  is  a  maximum  such  that  the 
contribution  to  the  pressure  field  of  the  "image"  is 
completely  suppressed,  then 


^(7) 


=  j  JJyr)r( 


{egn  2.41) 


and  the  magnitude  of  the  spectrum  is 


g(Y) 


(egn  2. 42) 


E.  THE  EFFECT  OF  ADDING  NOISE 

In  an  effort  to  simulate  a  realistic  ocean  environment, 
a  normally  distributed  narrowband  noise  field  based  on  a 
Rician  distribution  [Eef-  6:  p.  189]  was  adapted.  The 
resulting  acoustic  pressure  field  can  be  expressed  as 
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^(R)  =  £(R)  +  n 


(eqn  2.  43) 


where  jq  is  a  normally-distributed  random  function  with  a 
mean  cf  zero,  a  standard  deviation  of  one  and  possesses  both 
amplitude  and  phase.  The  noise  operates  independently  on 
both  the  amplitude  and  phase  components  of  the  pressure 
field.  Since  n  is  independent  of  range  and  wavenumber,  it 
is  treated  as  a  constant  by  the  Fourier  Transform. 

A  fourth  intention  of  this  thesis  was  to  discover  the 
degree  to  which  ncise  degrades  the  wavenumber  spectrum 
compared  with  the  corresponding  degradation  of  the  range 
domain.  To  produce  the  random  noise  fields  used  in  this 
research,  the  routine  listed  as  “dGNMl"  in  the  IMSL  library 
[Ref.  9]  was  called  upon  twice  to  generate  independent, 
pseudo-random  functions  which  interact  separately  and  simul¬ 
taneously  with  both  the  amplitude  and  phase  components  of 
?{R)  .  While  not  an  elaborate  scheme,  owing  to  time  and 
resource  constraints,  it  was  felt  that  this  modest  simu¬ 
lacrum  of  narrowband  noise  would  give  a  fairly  accurate, 
first  cut  "feel”  for  the  effects  of  noise  on  the  wavenumber 
spectrum.  The  reader  should  be  aware,  however,  that  no 
direct  consideration  cf  coherency  along  the  range  path  was 
taken  into  account  by  this  method. 
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Figure  2.2 


A  Classic  jP(R)|  vs. 


R  Curve 
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III.  EESOLTS  AND  CONCLO SIGNS 

A.  TBZ  FFT  ALGOBITHl 

The  research  for  this  thesis  was  conducted  entirely  with 
computer  algorithms  to  model  the  Lloyds  Mirror,  the  Fast 
Fourier  Transform,  the  range  windows,  waves,  and  the  narrow¬ 
band  noise.  This  was  necessary  since  the  Lloyds  Mirror  is 
cin  idealized  representation  of  a  situation  that  occurs  only 
rarely  in  nature,  and  at  that  is  limited  to  isospeed  surface 
ducts.  Also,  time  and  money  constraints  were  such  that  the 
use  of  computer  models  was  an  absolute  necessity.  For  these 
reasons  and  because  the  FFT  is  the  heart  of  the  Wavenumber 
Technique,  finding  a  reliable,  easy-to-use  computer  algo¬ 
rithm  to  perform  the  FFT  on  the  complex  pressure  function 
was  deemed  very  important. 

Initially,  the  IMSL  library  routines,  FFT2C  and  FFTCC 
[Ref.  9:  p.  232],  were  used  to  generate  the  pressure 
spectra.  Although  the  author  was  unable  to  ascertain  just 
how  Stamey  performed  the  FFT  on  his  data,  it  is  highly  prob¬ 
able  that  he  used  an  IMSL  library  routine.  While  in  nis 
thesis,  Blanchard  stated  that  he  used  the  FFTCC  routine. 
The  magnitude  of  a  typical  spectrum  generated  by  the  routine 
FFT2C  is  shown  in  Figure  3.  1. 

The  same  scenario  used  to  generate  Figure  3- 1  was  used 
to  produce  the  spectrum  shown  in  Figure  3.2,  but  the  routine 
FFTCC  was  utilized.  As  can  be  seen,  the  resulting  graphs  of 
the  pressure  spectrum  magnitude  are  practically  identical. 
Indeed,  the  set-up  of  the  pressure  fields  for  insertion  into 
each  routine  differ  cnly  in  that  FFT2C  requires  a  data  set 
consisting  of  an  integral  power  of  two  number  of  points, 
whereas  FFTCC  will  handle  any  number  of  points.  And,  wnen 
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the  Cooley-Tulcey  FFT  listed  in  Appendix  I  was  used  to  trans¬ 
form  the  same  pressure  field  as  was  utilized  to  produce 
Figures  3.1  and  3.2,  Figure  3.3  was  the  result.  This  algo¬ 
rithm  differs  somewhat  from  the  IKSL  routines  by  the  manner 
in  which  the  pressure  field  data  is  treated.  Both  FFT2C  and 
FFTCC  are  designed  to  handle  a  complex  array.  The 
Cooley-Tukey  FFT,  on  the  other  hand,  is  designed  to  trans¬ 
form  real  data;  but  the  pressure  field  is  complex.  To  get 
around  this  apparent  conundrum,  the  pressure  data  is  defined 
as  either  real  or  imaginary  by  its  placement  in  the  array. 
As  can  be  seen,  the  resultant  graph  of  the  spectrum  magni¬ 
tude  is  almost  indistinguishable  from  that  produced  by 
either  IMSL  routine-  In  conclusion,  it  does  not  matter 
which  FFT  routine  is  used  to  generate  the  pressure  spectrum, 
as  long  as  the  pressure  field  in  the  range  domain  is 
"fitted"  into  the  array  in  a  manner  suitable  to  the 
particular  algorithm  used. 


TABLE  1 

Critical  Values  Osed  in  the  Research 


K  (m“^) 

N 

r  (m) 

Cm 

3.0 

1024 

0.524 

1.172  X  10'^ 

2.0 

1024 

0.785 

7.816  X  10"^ 

1.0 

1024 

1.571 

3.906  X  10"^ 

It  is  important  to 

note  four  things 

about  the  FFT  at 

point. 

.  Since  the 

FFT  is 

given  only  a  f 

inite  sample  nf  a 

function 

whi  le 

the  theoretical 

Fourier  Transform 
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looks  at  the  complete  function,  the  computer  algo¬ 
rithm  is  designed  to  overcome  this  necessary  short¬ 
coming  by  assuming  the  submitted  finite  sample 
repeats  its  pattern  an  infinite  number  of  times. 
This  can  present  problems.  If  the  right  hand  side  of 
one  pattern  does  not  flow  smoothly  into  the  left  hand 
side  of  the  next  repetition,  as  is  shown  in  Figure 
3.4,  then  high  frequency  oscillations  known  as  the 
Gibbs  Phenomenon  can  be  introduced  into  the  spectrum. 
(See  Figure  3.5)  This  is  a  very  real  danger  where 
accustic  pressure  fields  are  concerned.  For  example: 
if  the  sample  size  does  not  include  the  entire 
pattern  (i.e.  ,  the  pressure  magnitude  has  not  been 
completely  attenuated  prior  to  the  endpoints  of  the 
sample) ,  this  unwanted  "jitter”  in  the  spectrum  can 
result.  To  avoid  this  situation,  not  only  was  a 
pressure  field  symmetric  about  tne  origin  (the 
source)  used,  but  also  a  Hanning  Window  was 
constructed  to  reduce  both  the  amplitude  and  the 
phase  components  of  the  pressure  field  to  zero  at  the 
endpoints,  thus  avoiding  the  Gibbs  Phenomenon. 
Figure  3. 6  shows  the  same  pressure  field  depicted  in 
Figure  3.4  combined  with  a  Hanning  Window.  Its  spec¬ 
trum  is  shown  in  Figure  3.7.  The  Gibbs  Phenomenon  is 
absent.  The  iorm  of  the  window  is  given  below.  The 
reader’s  attention  is  directed  to  Appendix  I  for  a 
view  of  the  manner  in  which  it  is  combined  with  the 
pressure  field. 

The  Hanning  Window; 


(eqn  3.1) 
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where  "NPTS"  is  the  number  of  data  points, 
and  I  =  1  ,2, 3, .. . ,NPTS. 

The  Hanning  Window  does  not  affect  the  acoustic  pres¬ 
sure  field  in  any  other  way;  the  interference  pattern 
in  the  range  domain  characteristic  of  the  Lloyds 
Kirror  for  an  Acoustic  Doublet  remains  the  same 
except  for  its  amplitude  near  the  end  points.  As 
shall  be  shown  later  in  this  chapter,  it  is  the 
pattern  of  the  entire  pressure  field,  not  simply  the 
amplitude  of  the  pattern,  on  which  the  P. T.  performs 
its  legerdemain. 

Because  the  vertical  wavenumber  beta  is  real  and  has 
physical  meaning  for  values  of  the  horizontal 
wavenumber  less  than  k,  only  those  values  of  the 
magnitude  of  the  acoustic  pressure  spectrum  —  cor re¬ 
sponding  to  <  k  were  retained. 

In  the  Wavenumber  Technique,  the  relationship  between 
the  horizontal  range  step  size  /^r  and  the  hori¬ 
zontal  wavenumber  step  size  is  extremely 
important.  Secalling  equation  2.36, 


Lyi^r  = 


(eqn  3.2) 


for  a  given  "N,”  the  term  on  the  right  hand  side, 
2  TT  /N  /  is  a  constant.  Consequently, 


Ar  «  (^y) 


(eqn  3.3) 


Simultaneously,  in  accordance  with  sampling  theory, 
in  order  to  prevent  aliasing  in  the  pressure  spec¬ 
trum,  there  is  an  upper  limit  on  the  size  of  /\r : 
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Ar  = 


{egn  3.4) 


^ r  is  chosen  such  that  it  is  much  less 
than  half  a  wavelength;  specifically,  <  X  /5  or 

X/6»  thus  a  good  sample  of  the  range  domain  is 
ensured.  However,  from  equation  2.36,  the  smaller 
Ar  becomes,  the  larger  A^  becomes  and  the 
coarser  the  sample  sizes  of  the  spectrum  become. 
Also,  to  ensure  adequate  samples  in  the  spectrum, 
examination  of  equations  2-25  and  2.36  and  [fief.  10] 


reveals  that 


Ar  = 


(eqn  3.5) 


where 


source  or  receiver. 


is  the  deeper  of  the  two  depths,  whether 
receiver.  So  there  is  a  delicate  balance 


to  maintain  between  the  sizes  of  Ar  and  At.  ir. 
the  spirit  of  compromise,  the  investigator  chose  to 
use  a  range  step  size  equal  to  a  quarter  of  a  wavel¬ 
ength  where 


X  =  — 

^  k 


(eqn  3.6) 


For  most  cases  investigated,  values  for  k  ,  N  ,  Ar 
and  At  are  summarized  in  Table  I.  As  can  be  seen, 
using  these  values  and  equation  3.5,  for  a  "Z"  value 
of  25  wavelengths,  the  maximum  size  At  can  assume 
and  still  ensure  an  adequate  sample  size  is 
At  '  0.0004  /m-  Clearly  the  sizes  of  At  shown 
in  Table  I  are  too  coarse  by  a  factor  of  approxi¬ 
mately  20.  Only  by  increasing  N  ,  the  number  of 
sample  points,  can  this  discrepancy  be  corrected. 
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However,  as  mentioned  earlier,  computer  resources 
were  limited.  The  effect  on  the  research  results  is 
not  significant;  nevertheless,  it  is  recommended  the 
reader  bear  this  limitation  rn  mind  during  the 
succeeding  pages. 

4.  Looking  again  at  Figures  3.1  through  3.3,  the  FFT  is 
the  curve  marked  by  the  solid  line,  while  the  theo¬ 
retical  Fourier  Transform  is  the  curve  marked  by  the 
dotted  line.  Notice  the  large  discrepancy  between 
the  theoretical  values  and  the  actual  values  in  the 
region  corresponding  to  small  vertical  wavenumbers. 
This  happens  to  correspond  to  a  region  of  the  spec¬ 
trum  characterized  by  more  rapid  fluctuations  in  the 
spectrum.  However,  the  size  of  beta  is  constant. 
Consequently,  fewer  samples  of  this  region  of  the 
spectrum  are  available  as  compared  with  that  region 
corresponding  to  high  beta  values.  Since  the  plot¬ 
ting  routine  utilized  is  interpolating  between 
sampled  points  of  the  spectrum,  the  graph  at  this  end’ 
is  less  smooth  and  results  in  a  marked  difference 
between  the  curves,  Where  the  curves  are  determined 
by  more  points  and  are  smoother,  the  discrepancy 
between  theory  and  computerized  reality  is  much  less 
noticeable. 

B.  SCDECB  DEPTH  DETEB3INATI0N 

Theory  predicts  that,  for  a  source  shallower  than  the 
receiver,  when  the  magnitude  of  the  pressure  spectrum  is 
plotted  as  a  function  of  the  vertical  wavenumber,  the  resul¬ 
tant  nodal  spacing  can  be  used  to  determine  source  depth  by 
either  consulting  a  chart  (see  Figure  3.8)  or  by  performing 
a  simple  calculation  (see  equation  2.31).  To  test  this 
prediction,  several  scenarios  were  designed.  Figures  3.9 


38 


REPRODUCED  AT  G  C  V  T  f’ rJ  r  i  r  t  J ';vp>rr4c 


through  3.11  are  typical  of  the  many  cases  run.  In  each, 
source  depth  is  less  than  receiver  depth.  For  each  case, 
careful  measurement  of  the  A/?  spacing  coupled  with  use  of 
equation  2.31  yielded  the  (already  known)  source  depth  with 
less  than  a  two  per  cent  error,  which  can  be  attributed  to 
human  error  in  the  measurement  of  A(0. 


TABLE  II 

Results  of  Source  Depth  Determination  Buns 


Figure  Nr. 

ZS  (m) 

ZR  (M) 

(Meas ' d ) 

Z(Calc ’d) 

3.9 

21.998 

47 . 124 

0.14  m’^ 

22.4  n 

3.10 

31.416 

47 . 124 

0.10  " 

31.4  '' 

3.11 

15.708 

47.124 

0.20  ” 

15.7  '■ 

3.12 

47 . 124 

21.998 

0.14  ” 

22.4  " 

3  . 13 

47 . 124 

31.416 

0.10  " 

31.4  " 

3  . 14 

47.124 

15.708 

0.20  ” 

15.7  " 

To  test  what  happens  when  the  receiver  depth  is  shal¬ 
lower  than  the  source  depth,  the  cases  depicted  in  Figures 
3.12  through  3.14  were  run.  The  scenario  is  the  same  as  tor 
Figures  3.9  through  3.11;  only  the  source  and  receiver 
depths  have  been  exchanged.  As  can  be  seen,  the  AA 
spacing  corresponds  with  the  receiver  depth.  This  is  again 
as  predicted  by  theory.  Table  II  summarizes  the  findings  of 
Figures  3.9  through  3.14. 
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The  scenarios  shown  in  Figures  3.9  through  3.14  all 
start  with  a  pressure  field  measured  horizontally  from  the 
source;  i.e.,  the  exact  range  to  the  source  is  known.  Tihat 
if  the  range  to  the  source  is  not  known,  or  only  minimal 
information  concerning  the  range  is  known?  In  this  case, 
only  a  portion  of  the  pressure  field  starting  at  seme 
initial  range  P.^,  can  be  sampled.  What  is  the  effect  on 
the  spectrum? 

To  answer  these  guestions,  a  range  window  was 
constructed  and  placed  at  varying  initial  ranges  from  the 
source.  As  Figure  2.2  illustrates,  as  range  from  the  source 
increases,  the  nodes  of  the  Lloyds  Mirror  interference 
pattern  in  the  range  domain  become  more  widely  spaced. 
Given  fewer  nodes  to  sample,  will  the  spectrum  be  the  same? 
In  accordance  with  theory,  the  spectrum  is  independent  of 
range  (see  equation  2.39)  and  so  should  be  the  same  wherever 
the  window  is  placed. 

Figures  3.15  through  3.17  illustrate  that  the  spectrum 
is  not  entirely  independent  of  range.  The  figures  utilize 
the  same  scenario  as  was  used  above,  differing  only  in  the 
range  from  the  source  at  which  the  sampling  of  the  pressure 
field  begins.  As  can  be  seen,  beginning  at  the  right  hand 
side  which  corresponds  to  high  vertical  wavenumber  values, 
or  the  low  horizontal  wavenumber  values  (see  equation  2.8)  , 
a  "washing  out"  of  the  spectrum  occurs,  increasing  as  the 
range  window  is  moved  further  from  the  source.  Again,  the 
theoretical  curve  is  marked  by  the  dotted  line. 

When  the  magnitude  of  the  pressure  spectrum  is  plotted 
as  a  function  of  the  vertical  wavenumber,  the  useful  portion 
of  the  resulting  graph  is  the  last  two-thirds  of  the  beta 
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the  following  criterion 
was  established:  the  nodal  spacing  of  the  spectral  plots 

was  no  longer  deter ninable  once  the  initial  range  was 

increased  to  approximately  three  times  the  source  depth.  In 
other  wordSy  minimum  knowledge  of  the  range  from  receiver  to 
source  must  be  available  to  ensure  adequate  sampling  of  the 
acoustic  pressure  field  in  the  range  domain.  This,  in  turn, 
will  produce  a  pressure  spectrum  of  such  a  quality  that 
source  depth  determination  can  be  made. 

TJhat  causes  this  not  entirely  unexpected  result?  As  the 
range  window  is  moved  further  from  the  source,  for 

qc  Zp  ,  those  samples  of  the  pressure  field  corresponding 
to  low  gamma  values  are  lost  to  the  spectrum  first.  low 

gamma  values  correspond  to  waves  arriving  at  high  angles 
with  respect  to  the  horizontal.  At  greater  receiver  ranges, 
?  more  closely  approximates  r  and  the  arrival  angles  become 
closer  tc  the  horizontal  (see  Figure  2.1). 

As  the  window  is  moved  further  from  the  source  and  fewer 
nodes  of  the  interference  pattern  are  encountered,  if  a 
wider  window  were  used,  could  this  ’’washing  out”  effect  in 
the  spectrum  be  minimized  or  even  eliminated?  In  the 
author’s  opinion,  it  would  be  minimized;  because  of  the 
reasoning  in  the  preceding  paragraph,  it  probably  would  not 
be  eliminated.  However,  a  wider  window  would  require  more 
sample  points;  because  of  equation  2.36,  one  cannot  merely 
make  [\z  larger.  Time  and  the  computer  resources  available 
to  this  investigator  were  limited,  and  the  test  cases  run 
were  already  using  the  maximum  available  array  sizes,  thus 
precluding  a  quantitative  illustralion  of  this  hypothesis. 
Succeeding  research  into  the  W.T.  might  include  ways  to  test 
this  point. 

In  conclusion,  the  Wavenumber  Technique  has  the  poten¬ 
tial  for  being  a  valuable  operational  tool  in  determining 
the  depth  of  an  acoustic  source  provided  the  receiver  is  at 
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a  greater  depth  than  the  source.  It  is  strongly  reccmmended 
that  further  research  into  use  of  the  H.T.  in  a  realistic 
ocean  environment  be  done. 

C.  THE  EFFECT  OF  SOEFICE  ROOGBNESS 

This  and  the  succeeding  section  were  purely  qualitative 
portions  of  the  research.  Therefore,  in  the  figures  the 
spectra  were  plotted  as  functions  of  the  horizontal 
wavenumber  for  ease  of  discussion  except  where  a  particular 
point  about  the  vertical  wavenumber  was  being  illustrated, 
and  the  theoretical  curve  was  omitted. 

As  described  in  Chapter  II,  the  effect  of  surface  rough¬ 
ness  is  to  suppress  the  contribution  or  the  image  to  the 
dipole  interference  pattern.  Moreover,  by  increasing  M,  the 
effect  is  to  suppress  the  image’s  contribution  to  the  pres¬ 
sure  field  in  the  range  domain  and  reduce  the  pressure  spec¬ 
trum  to  the  contribution  of  the  direct  path  wave  only  (see 
equation  2.41);  the  resultant  magnitude  of  the  spectrum  is 
inversely  proportional  to  the  vertical  wavenumber  (see  equa¬ 
tion  2.42).  As  Figures  3.13  through  3.21  illustrate,  this 
is  exactly  what  happens  in  the  FFT.  In  Figure  3.21,  an 
inverse  beta  curve  has  been  manually  superimposed  onto  the 
actual  curve  to  reflect  this  point. 

Looking  at  the  pressure  field  in  the  range  domain  ^(F) 
the  effect  of  surface  roughness  is  inversely  proportional  to 
the  wavelength  of  the  signal.  For  longer  wavelengths,  the 
Lloyds  Mirror  effect  is  more  tolerant  of  surface  roughness 
than  it  is  for  shorter  wavelengths.  Does  the  same  relation¬ 
ship  hold  for  the  pressure  spectrum?  Comparison  of  the 
respective  sea  states  in  Figures  3.22  and  3.23  with  Figures 
3.24  and  3.25  illustrate  that  it  does. 

This  investigator  also  looked  at  the  effect  of  varying 
the  range  window  in  the  presence  of  waves.  This  is  similar 
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to  what  was  looked  at  in  Section  B  of  this  chapter.  Results 
should  be  similar  and  for  the  same  reasons.  For  a  given  sea 
state,  source  and  receiver  geometry  and  wavelength,  the 
range  window  was  moved  successively  further  from  the  source 
(see  Figures  3.26  through  3.29).  As  can  be  seen,  the 
’’washing  out”  effect  cf  the  spectral  nodes  begins  to  effect 
the  lower  horizontal  wavenumber  values  first. 

D.  THE  EFFECT  OF  BOISE 

Random  noise  was  simulated  in  this  research.  Being 
independent  of  range  and  frequency,  it  was  treated  as  a 
constant  by  both  the  Fourier  Transform  and  the  FFT.  The 
investigator  wished  to  compare  qualitatively  the  effects  of 
equal  amounts  of  noise  on  the  pressure  field  in  the  range 
domain  and  in  the  wavenumber  domain.  To  do  this,  succes¬ 
sively  larger  values  of  the  ’’noise  factor"  were  used  in 
computing  ^  (see  equation  2.43).  Values  of  the  noise 
factor  used  in  the  research  were  limited  to  the  maximum 
amplitude  cf  the  pressure  field  ^(E) .  Since  dipole  radia¬ 
tion  is  predicated,  as  range  from  the  source  increases,  the 

maxima  of  the  interference  pattern  decrease  as  the  inverse 

_  2. 

square  of  the  range  (i.e.,  R  ).  Hence,  even  small  noise 
factors  can  have  a  devastating  effect  on  the  pressure  field 
as  the  receiver  is  stepped  out  in  range.  The  destructive 
effect  on  the  spectrum  of  successively  more  intense  amounts 
of  noise  at  varying  wavelengths  is  illustrated  in  Figures 
3.30  through  3.47 

Several  conclusions  concerning  the  pressure  spectra  can 
be  drawn  from  these  results: 

1.  For  equal  amounts  of  noise,  longer  wavelengths  are 
affected  less  than  shorter  wavelengths.  Compare 
Figures  3.31  and  3.32  with  Figures  3.43  and  3.44,  and 
Figures  3.34  and  3.35  with  Figures  3.46  and  3.47. 
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This  was  cot  a  surprising  result  since  the  sane  prin¬ 
ciple  holds  for  the  range  domain. 

2.  The  spectrum  is  affected  less  ty  noise  than  is  the 
range  domain.  Compare  Figures  3.30  through  3.35  and 
Figures  3.36  through  3.47  for  an  illustration  of  this 
point  at  two  different  wavelengths  and  various  noise 
levels.  This  was  a  rather  welcome  surprise. 

3.  If  the  magnitude  of  the  pressure  spectrum  is  plotted 
as  a  function  of  beta,  the  ^ /O  spacing  can  still 
be  detected  even  after  the  pressure  function  in  the 
range  domain  has  been  "swallowed  up”  by  the  noise 
field.  Compare  Figures  3.33  through  3.35  and  Figures 
3.42  through  3-44.  This  was  not  surprising  in  view 
of  2.  above. 

E.  SOHBATIOH 

All  of  the  foregoing  results  are  based  specifically  on 
the  lloyds  Mirror  for  an  Acoustic  Doublet.  This  is  a  highly 
idealized  and  artificial  environment.  Any  thought  of  imme¬ 
diately  applying  conclusions  reached  in  this  paper  to  a  real 
world  situation  is  premature.  However,  the  results  are- 
still  significant  if  only  for  supporting  the  statement  made 
by  each  preceding  investigator  that  the  Wavenumber  Technigue 
should  be  examined  very  closely  for  a  future  operational 
role,  especially  in  view  of  the  current  trends  in  source 
levels  and  ambient  noise  levels. 

Certain  fundamental  conclusions  regarding  the  author’s 
research  into  the  Wavenumber  Technigue  can  now  be  made 
within  the  confines  of  the  above  restriction: 

1.  Source  depth  can  be  determined  guickly  and  easily 
from  the  acoustic  pressure  spectrum  provided 
a)  the  receiver  is  deeper  than  the  source. 
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b)  some  knowledge  of  the  range  from  receiver  to 
source  is  available, and 

c)  the  magnitude  of  the  pressure  spectrum  is  plotted 
as  a  function  of  the  vertical  wavenumber. 

Source  depth  determination  in  even  a  noisy  environ¬ 
ment  is  possible.  While  the  introduction  of 
increasing  amounts  of  narrow  band  noise  adversely 
affects  both  the  range  domain  and  the  wavenumber 
domain,  the  pressure  spectrum  seems  able  to  withstand 
the  chaotic  destruction  longer  than  does  the  pressure 
field  in  the  range  domain. 

The  ability  to  determine  source  depth  is  adversely 
affected  by  the  height  of  surface  gravity  waves  since 
surface  roughness  reduces  the  Lloyds  Mirror  as  the 
sea  state  increases. 
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Figure  3. 1  Pressure  Spectrum  Osing  FFT2C 
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Figure  3.2  Pressure  Spectrua  Using  FFTCC 
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PRESSURE  SPECTRUM  VS  KZ 

SOURCE  DEPTH  =  20  m,  K  =  2.0/m, 
RECEIVER  DEPTH  =  50  m.  N  =  1024 
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Figure  3.3  Pressure  Spectrua  Osing  Cooley-Tukey  FFT 
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Figure  3.5  Pressure  Spectrum  Showing  Gibbs  Phenomenon 
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Figure  3.6 

Pressure  Field  Combined  Hith  a  Hanning  Window 
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RECEIVER  DEPTH  =  47.124  M,  K  =  2.0/M, 
N  =  1024 
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Figure  3.11  Graph  of  Pressure  Spectrum, 
Source  at  15.7  Meters 
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Figure  3-12  Graph  of  Pressure  Spectrum, 
Receiver  at  22-0  Meters 


PRESSURE  SPECTRUM  MRCNITUDE 


f?EPRODUCED  AT  G  O  V  T  t>  rg  M  r  t  J  T 


PRESSURE  SPECTRUri  VS  KZ 

SOURCE  DEPTH  =  47.124  M,  K  =  2.0/M, 
N  =  1024 


Figure  3.14  Graph  of  Pressure  Spectrua, 
ReceiTer  at  15.7  Meters 
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PRESSURE  SPECTRUM  VS  KZ 
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Figure  3. 15  Pressure  Spectrum, 
Range  Hindow  Set  at  47.1  fleters 
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Figure  3.18  Pressure  Spectruni,  Sea  State  0 
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Figure  3.19  Pressure  Spectrum,  Sea  State  2 
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Figure  3.20  Pressure  Spectrua,  Sea  State  3 
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SOURCE  DEPTH  =  21.Q<?8  M,  N  =  1024 
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Figure  3.24 
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SOURCE  DEPTH  =  21.998  M,  N  =  1024 
RECEIVER  DEPTH  =  47.124  M 


Figure  3.25  Pressure  Spectrum,  Sea  State  3,  K  =  2.0 
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Figure  3. 26  Pressure  Soectrum, 
SS  2/  Range  Window  Set  at  50  Meters 


Figure  3.28  Pressure  Spectrum, 

SS  2,  Range  Window  Set  at  200  Meters 
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Figure  3.29  Pressure  Spectrum, 

5S  2,  fiange  Hindov  Set  at  300  Meters 
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Figure  3.30  Pressure  Field,  K  =  1.0,  jj^  -  0.005 
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Figure  3.32  Pressure  Spectrum  vs.  Beta 
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Figure  3.34  Pressure  Spectrum  vs.  Gamma 
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Figure  3.35  Pressure  SpectruQ  vs 
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Figure  3.38  Pressure  Spectrum  vs.  Beta 
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Figure  3,40  Pressure  Spectrum  ts-  Gamma 
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Figure  3.41  Pressure  Spectrum  vs.  Beta 
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Figure  3.43  Pressure  Spectrum  vs.  Gamma 
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Figure  3.47  Pressure  Spectrum  vs 


APPENDIX  A:  LLOYDS  MIRROR  PRESSURE  FIELD  SOURCE  CODE 
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